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ABSTRACT 

Ultraluminous supersoft X-ray sources (ULSSS) exhibit supersoft spectra with 
blackbody temperatures of 50-100 eV and bolometric luminosities above 10^^ erg 
s~^, and are possibly intermediate mass black holes (IMBHs) of > IO^Mq or 
massive white dwarfs that are progenitors of type la supernovae. In this letter we 
report our optical studies of such a source in M81, M81-ULS1, with HST archive 
observations. M81-ULS1 is identified with a point-like object, the spectral energy 
distribution of which reveals a blue component in addition to the companion of an 
AGB star. The blue component is consistent with the power-law as expected from 
the geometrically-thin accretion disk around an IMBH accretor, but inconsistent 
with the power-law as expected from the X-ray irradiated flared accretion disk 
around a white dwarf accretor. This result is strong evidence that M81-ULS1 is 
an IMBH instead of a white dwarf. 

Subject headings: Galaxy: individual(M81) — X-rays: binaries 



1. INTRODUCTION 



Luminous supersoft X-ray sources were discovered with the Einstein Observatory and 
were established as an important new class of X-ray binaries based on ROSAT observations 
of 18 sources in the Milky Way and the Magellanic Clouds (Kahabka & van den Heuvel 
1997 and reference therein). These sources have extremely soft spectra with equivalent 
blackbody temperatures of 15-80 eV and are highly luminous with bolometric luminosities 
of 10^^ — 10^^ erg s~^. They are thought to be white dwarfs that are steadily or cyclically 
burning hydrogen accreted onto the surface. Remarkably, the accretion rates must be in a 
narrow range of 1 — 4 x 10~'^ MQyr~^ , and the resultant luminosities are below 10^^ erg s~^ 
as observed. White dwarfs with steady nuclear burning are promising Type la supernova 
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progenitors because, unlike explosive nova events, they can retain the accreted matter and 
their mass can increase until it approaches the Chandrasekhar limit. 

More and more supersoft sources have been discovered in distant galaxies with the 
advent of the Chandra and XMM-Newton X-ray observatories. These distant supersoft 
sources, as compared to the canonical ones, are generally hotter, more luminous and often 

associated with spiral arms, suggesting that some are young and massive systems (Di Stefano 
& Kong 2004). Some have bolometric luminosities, as derived from blackbody models, far 
above the Eddington luminosity {LEdd = 4:7rGMmpc/aT ~ 1.3 x lO'^^(M/A'/0) erg s^"*^) for 
white dwarfs, which we call ultraluminous supersoft sources (ULSSS), or ULS as a more 
compact name. One prototype ULSSS is MlOl-ULXl, which showed supersoft spectra of 
50 ~ 100 eV, 0.3-7 keV luminosities of 3 x 10^° erg s~^ and bolometric luminosities of ~ 10^^ 
erg s~^ during its 2004 outburst (Kong et al. 2004). While the supersoft spectrum can 
be explained by a white dwarf burning accreted materials on its surface, the white dwarf 
can not explain the luminosities that are extremely super- Eddington for white dwarfs of 
< IAMq. On the other hand, an intermediate mass black hole (IMBH) of > lO'^AfQ, as 
Kong et al. suggested, can naturally explain both the high bolometric luminosities and the 
supersoft spectrum given the scaling relation Tj„ oc M~^/^ between the black hole mass and 
the accretion disk inner edge temperature. 

Intermediate mass black holes (lO^-lO^M©) have played an important role in the hier- 
archical merging scenario of galaxy formation (Madau & Rees 2001). They have long been 
searched for in the cores of globular clusters, but iron-clad evidence is lacking (for a review, 
see Miller & Colbert 2004). Recently, much interest has been directed toward ultralumi- 
nous X-ray sources (ULXs), which have been found to exist outside the nuclear regions with 
Lx ~ 10^^-10^^ erg/sec in many nearby galaxies (Miller & Colbert 2004). Such luminosities 
require IMBHs of > 10^ M© if ULXs radiate at 10~^ Eddington level as do many active 
galactic nuclei and Galactic X-ray binaries. One type of evidence for IMBHs comes from 
the X-ray spectroscopy that reveals a cool accretion disk component with the inner edge 
temperature of a few xlO eV, which corresponds to an IMBH of a few xlO^ Mq if the inner 
edge of the accretion disk corresponds to the last stable orbit (Miller et al. 2004). However, 
in many cases the cool component is dominated by a hard power-law component. This hard 
component resembles those from the hot corona in stellar black hole X-ray binaries that will 
erode and envelope the inner part of the accretion disk, raising doubts against the IMBH 
interpretation and suggestions that ULXs may be black holes of 20-100 Mq in very high state 
(Done & Kubota 2006; Soria 2007). Indeed, ULXs can, from a theoretical perspective, be 
stellar mass black holes with beaming effects or radiating at super-Eddington levels with the 
photon-bubble instability operating in a radiation pressure dominated magnetized accretion 
disk (Begelman 2002). 
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Ultraluminous supersoft sources, with bolometric luminosities > 10^^ erg s~^ and su- 
pcrsoft spectra of 50 — 100 eV, are more promising IMBH candidates (Kong et al. 2004, 
2005). As compared to ULXs, ULSSSs exhibit spectra that do not show a hard power- law 
component, making it (more) reasonable to link the temperature to the inner edge of the 
accretion disk, thereby strongly suggesting an IMBH as the accretor. ULSSSs can be the 
ultraluminous version of the canonical supersoft sources as ULXs can be stellar black holes. 
However, given the nature and dimension of the surface nuclear burning for white dwarfs, it 
is unlikely to make apparent super-Eddington emitters with beaming effects or the photon- 
bubble instabihty in magnetized accretion disks as proposed for stellar black holes in ULXs. 
Of course, applying white dwarf atmosphere models to supersoft source spectra may reduce 
their bolometric luminosities by a factor of ten (van Teeseling et al. 1996), hence alleviating 
the super-Eddington problem for some ULSSSs, but it is unlikely to push the most luminous 
ones below the Eddington limit. 

Swartz et al. (2002) observed the nearby spiral galaxy M81 with Chandra ACIS, and 
discovered an intriguing ULSSS in the bulge at R.A. =09:55:42.2, Decl. =69:03:36. 5, which 
we designate as M81-ULS1 in this paper. Liu (2007) studied the spectral properties of this 
ULSSS with 17 Chandra ACIS observations, and found its spectra persistently supersoft 
in six years despite dramatic flux changes. The spectrum from the longest observation 
can be well fitted in 0.3-2 keV by an absorbed blackbody model for a nuclear burning 
white dwarf, with xl/dof = 1.196/38, uh = 8.6 ± 0.9 x 10^° cm-^, /cT = 73 ± 1.5 eV, 
Lx(0.3 - 2keV) = 3.2 x 10^^ erg s'^, and L^^i = 2.5 x 10^^ erg s'^ The spectrum can be 
equally well fitted by a multicolor disk model, with xl/dof = 1.073/38, uh = 10.5±0.9 x 10^° 
cm-2, kTin = 83 ± 1.4 eV and Lx(0.3 - 2keV) = 3.2 x lO^^ erg s'^. The inner disk radius 
derived from the model normalization is Rin ~ 24500/(cosi)^/^ km, which corresponds to an 
IMBH of ~ 2700/ {cosiy/^MQ, assuming corresponds to the last stable orbit. Thus, the 
X-ray spectra are suggestive of M81-ULS1 being either a white dwarf or an IMBH. 

Observations in other wavelengths can shed light on the nature of this class of X-ray 
sources. We have embarked on an effort to study such sources in nearby galaxies with HST 
observations. In this letter we report the optical studies of M81-ULS1 with HST archive data. 
§2 describes the HST archive data sets utilized, the optical identification and photometric 
measurements of ULSl. The photometric measurements of ULSl are interpreted in fight 
of the white dwarf model and the IMBH model in §3. We discuss our results and their 
impfications in §4. The distance of M81 is taken to be 3.63 Mpc (// = 27.8 mag; Preedman 
et al. 1994) in this work. 
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2. ANALYSIS OF HST DATA 

HST has observed the sky region around ULSl several times (Table 1). These include an 
observation on 1996-01-01 in fihers WFPC2/F547M (3x30s) and WFPC2/F656N (3x600s), 
an observation on 2003-09-18 in filters ACS/WFC F658N (700s) and ACS/WFC F814W 
(120s), an observation on 2004-09-14 in ACS/WFC F814W (1650s), and an observation on 
2006-03-27 in filters ACS/WFC F435W (465s) and ACS/WFC F606W (470s). All data sets 
were downloaded from MAST and calibrated on-the-fly with the best calibration files as of 
November 7, 2007. 

To identify the X-ray sources on the HST images, we register both the HST images and 
the Chandra images onto the 2MASS reference frame. This was achieved by identifying 11 
2MASS sources on the ACS/WFC F814W image taken in 2004, and two 2MASS sources on 
the combined Chandra image. After both were registered onto the 2MASS reference frame, 
the X-ray image was overlayed onto the ACS/WFC F814W image. Six X-ray sources were 
coincident with optical objects within error circles of 0'.'5, including ULSl that is identified 
with a point-like object. This object appeared in all HST observations (Figure 1), and is 
located in the smooth bulge of M81 with the closet object of equal brightness in F814W 10" 
away. 

The photometry for this point-like object was calculated with an aperture of 0'.'5 on 
WFPC2/PC1 and ACS/WFC images. The count rates and the corresponding magnitudes 
in the vegamag photometric system are listed in Table 1 for the counterpart in the seven 
observations. We calculate the magnitudes for two Ha filters with SYNPHOT assuming an 
Ha emission line atop the continuum of the Bruzual stellar spectrum bz_55; this spectrum, 
with an intrinsic color F547M — F814VF «i 1.3 as similar to observed, is normalized to give 
the observed F547M magnitude. The observations require an H^ emission line luminosity of 
(1.9 - 3.2) X 10^^ erg s"^ during the 1996 WFPC2/F656N observation, and (0.8 - 1.1) x lO^^ 
erg s~^ during the 2003 WFC/F658N observation. 



3. INTERPRETATION OF HST OBSERVATIONS 

The optical counterpart of ULSl, placed at the distance of M81, has an absolute mag- 
nitude of Mpsim ~ -7 mag, a color of F5A7M - FSMW ^ 1.3 as for K-M stars, and 
-F435W^ — F547M ^ —0.1 as for late B stars. The red color and absolute magnitudes in 
F547M and F814W might be obtained from a massive old globular cluster of ~ 10^ K-M 
dwarfs. However, the same globular cluster can not explain the late B-like blue color of 
F4:35W — F547M —0.1. Furthermore, the variabifities in the Ha and F814W bands 
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suggest that the counterpart is a single object. For single objects, such a brightness and 
red color in F547M and F814W arc not seen in dwarf /giant /supergiant stars as tabulated 
in Allen's Astrophysical Quantities (Cox 2000), but can be seen in asymptotic giant branch 
(AGB) stars that arc burning hydrogen/helium in shells and shine as a blackbody. We note, 
however, the object is considerably bluer than AGB stars with the same F814W luminosities, 
and there must be a bluer component in addition to the AGB star to account for the ob- 
served colors and magnitudes simultaneously. To quantify the AGB star properties in term 
of the core mass (Mc), we adopt the luminosity core-mass relation L — 59250(Mc — 0.522) 
(Paczynski 1970), the core-mass radius relation R/Rq = RqM^-^ /{l+4:M^)+0.5 {Rq = 4950; 
Rappaport et al. 1995) and Stefan-Boltzman's law L — AnR'^aT'^. 

A natural explanation for the blue component is that it originates from the accretion 
disk. If the accretor is a white dwarf, the accretion disk is illuminated by the soft X- 
ray emission from the white dwarf surface. Given the high X-ray luminosity, such a disk 
shall flare up geometrically, intercept a large fraction of the X-ray emission and re-emit 
in lower energies (Popham & Di Stefano 1996). The emergent spectrum will take a form 
of in the optical (FKR 2002). We calculate the F547M, F435W, 606W, and F814W 
(2004) magnitudes with SYNPHOT for an AGB star plus a powcrlaw of z/~^, for different 
core masses and relative contributions from the two components, and compute the quantity 
Am^ = S(mobs — rrimodeiY as a sum of the squared differences between the observed and 
model magnitudes. By minimizing Am^, we found at Am^ = 0.4 the best combination of 
the powerlaw component and an AGB star, which has Mc — 0.77Mq, R — 639i?0 and 
T — This model, as plotted in Figure 2, can not explain the observations because 

it does not give enough blue light in F435W. The blackbody spectral fit to the X-ray data 
(Liu 2007) is overplotted for comparison; clearly, a white dwarf itself would be too dim in 
the optical to explain the observations. 

If the accretor is an IMBH and thus the high X-ray luminosity is less than a few per- 
cents of the Eddington luminosity, the accretion disc will be a standard geometrically-thin 
optically-thick one (Shafee et al. 2007). The spectrum of such a thermal disk is convention- 
ally described by a multi-color disc (MCD) model (Mitsuda et al. 1984). The power-law 
component with F^, oc u^^^ for the MCD spectrum may extend down to the optical if the 
outer disk radius is large. For ULSl, this means Rgut ^ which is reasonably small 

(~ 4O/(cosi)^/^i?0) as compared to the giant star size or the Roche lobe size of the black 
hole. With similar procedures as in the white dwarf model, the four observed magnitudes 
determine a best combination, at Arm? — 0.16, of the v^/^ powerlaw component and an AGB 
star, which has Af, = IMMq, R = ISOli?© and T = 2277X. The model, plotted in Figure 3 
with the MCD model fit (Liu 2007) overplotted for comparison, can explain all four observed 
magnitudes. Remarkably, the required u^^^ power-law component in the optical is consistent 
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to within 30% with the MCD fit to the X-ray spectrum. 

The properties of the companion AGB star may differ considering the uncertainties in 
the theory and the observations. In particular, the analytical form for the core-mass radius 
relation for AGB stars, adopted for the convenience of calculations, tends to overestimate 
the radius (Rappaport et al. 1995). To test how this affects our results, we re-calculate 
the models by setting Rq = 3300, i.e., reducing the radius by 33%. Under this simplistic 
modification, we find an AGB star with = 0.60Mq, R = 2197?e and T = 3119fs: for 
the white dwarf model, and an AGB star with = 0.93Mq, R = 595Rq and T = 2859/^ 
for the IMBH model. The white dwarf model has a Am^ = 0.42 and cannot explain all the 
observations because it docs not give enough blue light. The IMBH model can fit all the 
observations simultaneously, but to a less satisfactory level (Am^ = 0.26) as compared to 
the model with i?o = 4950 (Am^ = 0.16). The newly computed AGB star is hotter and 
more compact. It has a less massive core (O.93M0), which is consistent with the final core 
masses of AGB stars with heavy stellar wind (Weidemann 2000). If the AGB star is at its 
final stage, its core mass indicates an initial mass greater than 5Mq and a lifetime shorter 
than 0.1 Gyrs (Vassihadis & Wood 1993). 



4. DISCUSSION 

Ultraluminous supersoft sources (ULSSSs), as X-ray spectroscopy suggests, are possibly 
the long-sought intermediate mass black holes (IMBHs), or massive white dwarfs that are 
progenitors of type la supernovae. We note that the accretion disk behaves quite differently 
in presence of an IMBH accretor or a white dwarf accretor for the same X-ray luminosities, 
and observations in optical and infrared arc capable to pick up the differences. In this 
Letter, we report the optical studies with HST archive observations of M81-ULS1, a ULSSS 
in the M81 bulge. To summarize, the spectral energy distribution of M81-ULS1 in the optical 
reveals a blue component in addition to the companion of an AGB star. The blue component 
is consistent with a, oc v^l^ power-law as expected from the geometrically-thin accretion 
disk around an IMBH accretor, but inconsistent with a F,^ oc pow-law as expected from 
the X-ray irradiated flared acrretion disk around a white dwarf accretor. This is strong 
observational evidence that M81-ULS1 is an IMBH instead of a white dwarf. 

The HST archive observations, however, were taken at different times spanning six years. 
This may affect our conclusion because the differences between different bands could result 
(partly) from the intrinsic variabilities of the accretion disk and the AGB companion. The 
accretion disk may have changed as suggested by the X-ray variability over the years, while 
AGB stars are known to be long period variables with amplitudes of up to AV^ ~ 2.5. The 
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marginally significant change in the F814W magnitudes from two observations separated by 
one year may have reflected such variabihties. Simultaneous observations in multiple bands 
are thus required to make clear whether and how much the spectral energy distribution in 
this work was contaminated by the variabilities. While such simultaneous observations may 
change the details of the AGB companion and the accretion disk, they probably will not 
change the main conclusion that the blue component comes from a oc v^l"^ powerlaw 
instead of a Fj, oc power. This is because the distinguishing power of the optical data 
mainly comes from the blue color between F435W and F606W (or F547M) as illustrated in 
Figures 2 and 3, yet F435W and F606W were already observed simultaneously. 

The conclusion would be invalidated if the optical counterpart were an interloper close 
to M81-ULS1 by projection. Such projection, however, is unlikely because the probability 
of finding such an AGB star in a (X'S error circle by chance is quite small (< 0.003), as 
calculated from the the error circle size and the fact that no stars of equal brightness were 
found within 10'.'. Also, the prominent emission is suggestive of the photoionization of 
surrounding materials by the soft X-ray emission from ULSl, and the Ha variation between 
two observations is somehow expected due to the long term X-ray variability of ULSl (Liu 
2007). Furthermore, a companion of an AGB star with a supergiant size, if overflowing or 
underfilling its Roche lobe, implies a period of 30 years or longer, which is consistent with 
the failure to detect any period from the six year Chandra observations (Liu 2007). 

We would like to thank Douglas Swartz, Scott Kenyan and Jeff McClintock for helpful 
discussions. JFL acknowledges the support for this work provided by NASA through the 
Chandra Fellowship Program, grant PF6- 70043. RD acknowledges the supports from grants 
NAG5-10705, AR-10948.019 and NNX07AH64G-R. 
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Fig. 1.— The optical counterpart for M81-ULS1 on HST images in F547M (1996), F656N 
(1996), F658N (2003), F814W (2004), F435W (2006) and F606W (2006). The WFPC2 
F547M/F656N and ACS/WFC F435W images are plagued with cosmic rays. All images are 
aligned with the counterpart at the image centers. The error circles have a radius of O'.'S. 
North is up in all images. The counterpart also appeared on the ACS/WFC F814W (2003) 
image, which was very noisy due to the short exposure and not shown here. 



- 10 - 



10" 



o 

B 

o 

\ 

M 

u 



10" 



10" 



10" 




0.1 



1 10 

photon energy (eV) 



100 



1000 



Fig. 2. — Fit for the optical observations with an AGB star plus an accretion disk illuminated 
by X-ray emission from the surface of a nuclear burning white dwarf. The X-ray illuminated 
accretion disk radiates as a power of F^, oc z/~^. The unabsorbed and absorbed spectral fits 
as a 73 eV blackbody for the white dwarf X-ray emission are plotted for comparison. 
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Fig. 3. — Fit for the optical observations with an AGB plus the power-law component of 
a multicolor disk around an IMBH. The power-law component has a form of oc u^^^. 
The multicolor disk has an inner edge temperature of 83 eV as from X-ray analysis, and its 
unabsorbed and absorbed spectra are plotted for comparison. 
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Table 1. HST observations for M81-ULS1 



ID 


Filter 


ExpT 


DATE 


Count Rate 


VEGAMAG 


U32L0102/5/8 


WFPC2/F547M 


3x30 


1996-01-01 


0.67 ±0.17 


22.1 ±0.2 


U32L0101/4/7 


WFPC2/F656N 


3x600 


1996-01-01 


0.09 ±0.01 


20.2 ±0.1 


J8MX18010 


WFC/F658N 


700 


2003-09-18 


8.7 ±0.5 


20.0 ±0.05 


J8MX18AKQ 


WFC/F814W 


120 


2003-09-18 


82.4 ±5 


20.7 ±0.06 


J90LA5010 


WFC/F814W 


1650 


2004-09-14 


73.3 ± 3 


20.8 ±0.04 


J9EL28A9Q 


WFC/F435W 


465 


2006-03-27 


29.9 ±0.5 


22.09 ±0.02 


J9EL28AAQ 


WFC/F606W 


470 


2006-03-27 


58.7 ±1.7 


21.98 ±0.03 



Note. — The columns are (1) exposure ID, (2) filter, (3) total exposure in seconds, (4) 
observation date, (5) count rate in e~/s for the counterpart, and (6) magnitudes in the 
VEGAMAG photometric system. 



